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Experimental evidence from both spin-valve and quantum transport measurements points towards 
unexpectedly fast spin relaxation in graphene. We report magnetotransport studies of epitaxial 
graphene on SiC in a vector magnetic field showing that spin relaxation, detected using weak- 
localisation analysis, is suppressed by an in-plane magnetic field, B\\ , and thereby proving that it is 
caused at least in part by spinful scatterers. A non-monotonic dependence of effective decoherence 
rate on B\\ reveals the intricate role of scatterers’ spin dynamics in forming the interference correction 
to conductivity, an effect that has gone unnoticed in earlier weak localisation studies. 


Loss of quantum information carried in the phase and 
spin of electrons propagating in a disordered conductor is 
associated with decoherence that suppresses interference 
corrections to conductivity The fundamental re¬ 

lation between spin relaxation and the low-temperature 
magnetoresistance (MR) has been investigated theoreti¬ 
cally and experimentally in numerous disordered metallic 
and semiconductor structures. 

Recent studies of spin-valve w and quantum trans¬ 
port [Ml effects in graphene-based devices have re¬ 
turned unexpectedly fast spin relaxation in graphene, 
despite this one-carbon-atom-thin material bein g b illed 
as an ideal medium for spintronics applications [l3l - ll5] . 
Spin relaxation can be induced by either spin-orbit cou¬ 
pling or spin-flip scattering due to magnetic impuri¬ 
ties. In quantum transport, spin-orbit coupling inverts 
the conventional negative MR around zero perpendicular 
held, B±, (weak localization, WL) to positive MR known 
as weak anti-localisation, but this effect has never been 
observed in pristine exfoliated or epitaxial graphenes 
oEm. By contrast, spin-hip scattering on magnetic 
impurities washes out quantum interference effects [il- 


1^ . including weak (anti-)localization. Spin-flip scatter- 
ing in disordered metals has been shown 2l-25| to raise 
the electron decoherence rate t~^{T) above the value ex¬ 
pected from inelastic scattering on thermal charge huc- 
tuations 12611 . 


T-^ = {kBT/h){peyh) ln(/r/2e>) = AT, (1) 


where p is resistivity. Residual decoherence in the limit 
T —>■ 0 is then determined only by the electron spin re¬ 
laxation time, Ts- Spin-hip decoherence was previously 
observed in exfoliated graphene nn. Applying a large 
in-plane magnetic held, Ry > Bp = ksT/gips, polarises 
impurity spins with g-factor gi, suppressing spin-hip scat¬ 
tering and prolonging phase coherence 0. For magnetic 
helds too small to polarize impurity spins, iJ|| < Bp, one 
might expect minimal effect on decoherence. 

In this paper we show both experimentally and theo¬ 
retically that such small in-plane helds do in fact have 
a discernable effect on phase coherence, but the effect 
is opposite to that observed for larger iJ||. The mea¬ 
surement is performed on epitaxial graphene grown on 
silicon carbide (SiC/G), using curvature of the B± MR 
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peak to quantify the electron decoherence rate. Apply¬ 
ing an in-plane magnetic field first broadens the MR 
peak slightly (enhances decoherence), before the sharp¬ 
ening effect due to impurity polarization sets in. This 
magnetic field dependence shows that the observed de¬ 
coherence is caused by spin-flip scattering rather than 
other dynamical sources of decoherence, such as external 
noise due to external two-level systems [l^. The non¬ 
monotonic dependence of decoherence rate on has not, 
to our knowledge, been discussed in previous work. It 
is a generic feature of quantum transport in disordered 
conductors that can be attributed to the precession of 
impurity spins at the frequency difefrent from the spin 
precession of mobile electrons. 

When electron (e) and impurity (i) g-factors differ, 
the difference between their spin precession frequencies, 
<^e|i = leads to a random variation of the 

spin-dependent scattering conditions for electron waves 
retracing the same closed diffusive trajectory in clock¬ 
wise and anti-clockwise directions, whose interference 
forms the quantum correction to conductivity. The non¬ 
monotonicity is characterized by a magnetic field scale, 
= hT~^/\ge — gi\g-B, above which the decohering effect 
of the in-plane field (separating precession frequencies for 
impurities and conduction electrons) is overcome by po¬ 
larization of impurity spins. The g-factor of magnetic 
scatterers can thus be determined by fitting the temper¬ 
ature and i?|| dependence of MR curvature to theory de¬ 
veloped below. 


Magnetotransport measurements were performed on 
a large area (160 /im x 35 ^m. Fig. [Ha)) n-doped 
SiC/G Hall bar, encapsulated in a polymer to improve 
temporal stability and doping homogeneity and exposed 
to deep-UV light to reduce carrier concentration (n = 
10 ± 1 X 10^^ cm“^) [ 2 ^. The device was measured in 
a dilution refrigerator equipped with a two-axis magnet, 
allowing independent control of B± and Hy llj. Aver¬ 
age MR measurements are not obscured by mesoscopic 
conductance fluctuations due to the large sample size. 


Figure (Hb) shows the characteristic negative MR of 
weak localisation, measured for i?|| = 0. As expected, 
MR is sharpest at the lowest temperatures where ther¬ 
mal charge fluctuations are minimized [Eq. (|H)]. Rais¬ 
ing the in-plane field to H|| = 1 T yields a significantly 
sharper (and higher) MR peak due to the partial polar¬ 
ization of impurity spins at this field [Fig. (He)]. Even in 
the raw data of Fig.[Tl(c), however, the non-monotonicity 
that is the primary subject of this paper can be clearly 
seen: the MR peak for Hy = 0.3 T is broader than the 
trace at Hy = 0, despite the small but non-zero impurity 
polarization at this low in-plane field. 


To make further progress, the decoherence time. Tip, 
is quantified using the expression for the curvature k of 


magneto-conductivity around B± = 0, 
d'^a 


dBl 


167r 

(DTp\ 

B_l=0 ^ ^ 
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( 2 ) 


that comes from the basic functional form of WL [ll- 
Q. Curvature is extracted from a parabolic fit to the 
average of 10 measurements of p{B±) covering the range 
of < 0.5 mT. The resulting temperature dependence 
of t~^{T) at Hy = 0 clearly shows the linear scaling 
expected from Eq. O [Fig.ura)]. 

The slope A of the temperature dependence in Eq. o, 
estimated using p^x = 1500 ± 35 H which we deter¬ 
mine from the measured resistance and the sample as¬ 
pect ratio yields A^. = 16.4 K“^ns“^ as compared to 
A Ri 31 K“^ns“^ fitted to experimentally-determined 
[Fig- EDS')]- This difference can be reconciled 
in two ways: (i) To note that the large-area epitax¬ 
ial graphene monolayers usually contain bilayer inclu¬ 
sions, which we also identified in the device used here, 
[Fig. (Ha)], and which have much higher conductivity 
than that of the monolayer material mil , hence, re¬ 
ducing the effective length of the Hall bar. (ii) To treat 
parameter A as an empirical factor 0. For the purposes 
of this paper we simply rescale the effective sample area 
to force the slope of t~^{T) to match that predicted by 
Eq. (|T|), giving pxx ~ 2800 H, mean free path ^ 26 nm, 
diffusion constant D = 131±10 cm^/s, and the slope 
As Ki 22 K“^ns“^. The black dashed line in Fig. |D[a) 


(a) 



(b) (C) 



FIG. 1. (a) Layout of the sample and an optical micrograph 
[ 3^1 showing predominantly monolayer graphene with bilayer 
inclusions; (b) Magnetoresistance sharpens with decreasing 
temperature (B|| = 0); (c) The sharpness of the magnetore¬ 
sistance curve is non-monotonic in i3|| (T = 450 mK). 
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(a) 


(b) 




FIG. 2. (a) Dephasing rate as a function of temperature at 
By = 0 extracted from the raw (blue circles) and rescaled 
(red circles) magnetoresistance. The blue dotted line is the 
best linear fit to the raw data. The slopes of the dashed lines 
are given by Eq. 0, and the offset constants are chosen 
arbitrarily for ease of comparison with the data. The solid 
red line is the theory developed in this work, (b) Dephasing 
rate extracted from the rescaled magnetoresistance plotted as 
a function of By at T = 450 mK. Roman numerals denote 
different scattering regimes. The dotted line is a guide to the 
eye. 


shows the plot of eq. m for the rescaled geometry, which 
closely matches the rescaled data and extrapolates to a 
finite value corresponding to the decoherence time ^ 20 
ps as r —>■ 0. Data presented in the remainder of this 
paper is based on the rescaled area, but we point out 
that the alternative approach (ii) to treat the parameter 
A does not alter the qualitative result of the following 
analysis of the quantum transport data. 

Figure [UJb) shows extracted from the curvature of 

rescaled magneto-conductance, for different values of By 
(T = 450 mK). The observed dependence t“^(B||) fea¬ 
tures three characteristic regimes. For the intermediate 
fields (II), the polarisation of spin-scattering impurities 
at ffi/isBii > fcsT, suppresses spin-flip scattering and 
decreases r“^. This prolongation of phase coherence by 
the in-plane field is a smoking gun for spin-flip scattering 
in the system. The suppression of is by a factor of 
approximately two before it begins to saturate at high 
field, indicating the entry into region III. The high-held 
saturation may be explained by hexural deformation of 
graphene out of the plane, resulting in randomly-varying 
flux from the in-plane magnetic held [s^ ; by g ~ 0 mag¬ 
netic moments that are not polarised by magnetic held; 
or by spin-orbit interaction. This regime will be the sub¬ 
ject of future work. 

In the analysis that follows, we focus on the lowest- 
held regime (I), where the decoherence rate extracted 
using Eq. ([2) shows a non-monotonic behaviour. Early 
studies of WL in dilute magnetic alloys demonstrated a 
‘spin-memory’ effect [ 1 , in the way the spinful scat- 
terers reduce the size of the WL correction to conduc¬ 
tivity. In summary, the exchange interaction between a 
spinful impurity and the electron spin s, via a perturba¬ 


tion Js ■ S(5(r — Ti), affects the electrons scattering in two 
ways: (a) by scattering an electron without hipping its 
spin and (b) by scattering with a spin hip. Although the 
spin hip process (b) always leads to the decoherence of 
electrons, process (a) leads only to additional scattering 
phases for a hxed spin conhguration of impurities (rela¬ 
tive to the electron spin). These phases would in general 
be the same for the two reversed sequences of visited scat- 
terers, hence, they would not suppress the interference 
correction to conductivity. However, if the spin conhgu¬ 
ration of impurities is randomly changing in time, e.g., 
by Korringa relaxation, electron waves travelling along 
closed paths in opposite directions and therefore arriving 
at the same scatterer at different times would experience 
randomly different conditions for a spin-conserving scat¬ 
tering (‘ergodic regime’), hence acquiring random relative 
|)lMse shifts that suppress their interference contribution 

In contrast to Korringa relaxation, spin precession 
of impurity spins does not necessarily randomise spin- 
dependent scattering conditions for diffusing electron 
waves. If electrons and scatterers have equal g-factors 
Si = He ~ 2 , their precessions in an external magnetic 
held are in phase and, because scattering conditions de¬ 
pend only on the relative orientation of the two spins, 
the ‘spin-memory’ regime would also persist in a hnite 
B||. For scatterers with gi ^ g^, on the other hand, spin- 
dependent elastic scattering amplitudes for electrons fol¬ 
lowing clockwise and anti-clockwise trajectories can be 
partly destroyed by even a small By. As sketched in 
Fig. [31(a), relative orientations of the spins of a scatterer 
and of the electron waves arriving along clockwise and 
anti-clockwise trajectories responsible for the interference 
correction to conductivity would typically deviate by an 
angle {g,, — gi)gBB\\T^p/h and become randomised when 
( 5 e — > hT~^. This effect should be strongest 

in systems where the difference between the electron and 
scatterer g-factors is largest (in particular, in the systems 
where these two have opposite signs). This leaves a range 
of magnetic held where precessional decoherence can fully 
develop (the ergodic regime) before the scatterers’ spins 
become closely aligned with the external held (By > Bt), 
thus diminishing the role of precessional dynamics. 


In WL theory, the difference between these ‘spin- 
memory’ and ‘ergodic’ regimes is accounted for by numer¬ 
ical factors that appear in the relations between dephas¬ 
ing rates of two-electron correlators, Cooperons, 

and the electron scattering rate from the spinful impu¬ 
rity, T-'^ = 2^^-fn^J'^S{S+l). Here, 7 = stands 

for the density of states of carriers and v is Dirac veloc¬ 
ity of electrons in graphene. The Cooperons are classihed 
according to their total spin S' = 0,1 and its projection 
M onto the direction of external magnetic held. These 
four Cooperons Cs,m are combined in the expression Q 
for the WL correction to conductivity 37| , 
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[^0.0 - C'l.o - Ci.i - Ci,_i]. 

Zirti 

In the spin-memory regime, Tg”Q = 2t“^ -1- ^ and 

, whereas in the ergodic regime [i^ all 
Cooperons decay with the same rate, Tg\j = T, 1 
Taking into account a possible difference between spin 
precession of the electrons and scatterers (which we do 
by changing the spin coordinates into the frame rotating 
with the frequency uje around the direction of external 
magnetic field Bi i), we find using diagrammatic pertur¬ 
bation theory [Sg that the two Cooperons Ci^±i with 
S' = 1 and M = ±1 are decoupled from all others and 
decay at the rate. 



{Sl)±{S,){l-2np{e^)) 
S(S + 1) 


T ’ 


where £ip = eigifisB/2 takes into account that spin-flip 
scattering requires energy transfer to the scatterer. At 
the same time, spin precession mixes the two Cooperons, 
Co,o a-nd Ci^o- This mixing generates combined modes, 
which relax with the decoherence rates 


'o.± 


(S^) - (S^)[nF(£+) - nF(£_)] 1 1 

S{S+1) \ " 


- {9e - . 


Note that, at S = 0, {Sz) = 0 and (S^) = \S{S -I-1), so 
that = '^oT- “ %+ = ^ 

corresponding to the relaxation rates of triplet and sin¬ 
glet Coperons in the spin-memory regime. At i3|| > 

reflecting the restoration of phase coherence 
of diffusion of spin-polarised electrons. 


Together, the four Cooperon modes [37| yield: 


a=± 

-1 


- 


In -I- Aa In 

'^iv ^ 0 ,+ 

' 0 ,- 1 - ~ '^ 0 ,- 




and magneto-conductance (measured as a function of B± 
for fixed i3||). 


cr(B±,Bii) = cr(0,B||) -^JdeJ2 


x{Fi4^) + 
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, . ,11, S/4e 

F{z) = \i\z + ili{-F -)■, = — -, 
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where n'p{e) = dnF(s)de is a derivative of the Fermi 
distribution function, ip is the digamma function, 
is the momentum relaxation rate that determines the 


classical Drude conductivity and the diffusion coefficient, 
D = (note that for monolayer graphene studied in 

this work, D RilSl cm^/s), and is the intervalley scat¬ 
tering time [H, [s^] . 

A non-monotonic dependence of decoherence rate on 
the in-plane magnetic field follows from the calculated 
magneto-conductance curvature. 


Stt/i 


YX {-n'pisaXx 

a—± 


X 


/ Dri^a 

V ^/e 


2 

+ (^lT±l -'^T ^ )(^0,-S+'ro.-) 


Dto,+ 

hje 


Dtq- 

hje 


(3) 


Here, we used the expansion F{z) k, z^/24 -|- 0[zX- 


(a) 


(b) 



FIG. 3. (a) Illustration of the influence of precession of the 
impurity spin on the weak localisation effect, (b) Comparison 
between the experimental values of the curvature of magne¬ 
toconductance at T = 0.45 — 2.1 K with the theoretical val¬ 
ues calculated using the parameters obtained using a fit at 
r = 2.1 K. 

Eq. ([3|) gives a semi-quantitative description of the ob¬ 
served magneto-conductance curvature that captures its 
qualitative features over a wide range of temperatures 
0.45 — 2.1 K and in-plane magnetic fields (0-3 T) [Fig. 3b]. 
Here, the following protocol was used for a single fit that 
yielded all five curves shown in Fig. 3b: 

1. The temperature dependence of magneto¬ 
conductance curvature k at B\\ = 0 is used 

to extract ss 77 ns“^ and the parameter 
At ~ 24 ns“^K“^ from the fit to our theory to 
be compared with As ~ 22 ns“^K“^ from eq. 
O- The red line in Fig. [2^ is the theoretical 
temperature dependence, which nearly coincides 
with the straight dashed line given by eq. m and 
perfectly matches the experimental points. 

2. The dependence of k on B\\ at high temperature, 
r = 2.1 K [bottom dataset in Fig. 3(b)], is used 
to determine the average g-factor of the scatterers, 
which returned the value gt = —0.84 ± 0.25. 
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The parameters determined above are sufficient to cal¬ 
culate the i?||-dependence of the magneto-conductance 
curvature for lower temperatures, down to T =450 mK, 
which are then compared with the experimentally mea¬ 
sured values [Fig. 3(b)]. The theoretically calculated 
curvature (lines) captures the observed non-monotonic 
behaviour and a shift of the low-field anomaly towards 
smaller at lower temperatures [s^ . 

To conclude, we have demonstrated that the excess de¬ 
coherence rate, observed earlier at low temperatures in 
epitaxial graphene sublimated on SiC [ 1 , @ , is caused by 
the spin-flip scattering of electrons from spinful impuri¬ 
ties. These spinful scatterers have an average ^-factor 
Qi = —0.84 ±0.25, very different from the free-electron 
g-factor in graphene, ( 7 e ~ 2, which enabled us to observe 
the influence of precession of the impurities’ spins on the 
WL effect. The large (and negative) g-factor for the im¬ 
purity spin implies a strong atomic spin-orbit coupling 
in the magnetic moment formation. This may indicate 
that these spins are in the surface states underneath the 
graphene layer, possibly originating in Si substitutions of 
carbon atoms in the interfacial layer (^ . The presence 
of such spinful scatterers on the SiC surface, directly ac¬ 
cessible for graphene electrons, explains the short spin 
coherence length observed in spin-valve experiments on 
a similar material Q. 
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THEORETICAL MODEL FOR THE 
DECOHERENCE ANOMALY IN 2D MATERIALS 
WITH MAGNETIC DEFECTS IN THE IN-PLANE 
MAGNETIC FIELD 


To explain the theoretical basis for the anomaly in 
k(-B||) discussed in the main text of the Letter, we con¬ 
sider an electron wave propagating along a closed-loop 
trajectory, scattering from disorder, 

^(r) = ^C/J(r-rc/y) JSj •S(5(r-rj,j), (1) 

r)7,j rj-j 

Zeeman splitting for electrons and impurities are taken 
into account by 


( 1 ) 



C , ,(£,q\Q)',co) 
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where z-projection cr of electron’s spin takes values 
+V 2 (t or -I-) or —Y 2 (i or —); g-factor of the electron 
ge K, —2.002, and g-s = \e\h/2mc is Bohr magneton. 

In the WL theory [ij, the interference correction, 5a, 
to conductivity is expressed as in term of two-electrons 
propagators, ’Cooperons’ as 


= — [Co.o - Ci,o - ; 

Cs,M = 


C =-^VpT^ / de(rc{np{erT')C{e,q). (2) 


corresponding to the diagramatic perturbation theory 
represenation shown in Fig. 1. Here, Clebsch-Gordan co¬ 
efficients = (S', Af IY 2 , ct;V 2 , cr') select from the 

Cooperon matrix C = {Ca-a-';riri') the singlet (S = 0) and 
triplet (S = 1, M = —1,0,1) components defined in 
terms of the total spin carried by the two-electron prop¬ 
agator and its projection onto the external magnetic field 
B|j. Also, Ea = e — agigsB^ and n'p{e) = dnp{e)/de 
is a derivative of the Fermi distribution function. 

To mention, in graphene, the weak localisation be¬ 
haviour may be additionally affected by the Berry phase 
of electrons arising from the sublattice composition of 
their wave functions on the honeycomb lattice of carbons 
0. This would lead to the antilocalisation behavior in 
structures with a long intervalley scattering time, turned 
into WL behaviour for long trajectories with the time of 
diffusive flight longer than the intervalley scattering time. 
In the materials studies here, intervalley scattering time 


FIG. 1. Line (1): The diagram for Kubo susceptibility cor¬ 
responding to the WL correction to the conductivity ex¬ 
pressed in terms of the Cooperon 

are current operators, the Keldysh components as well as spin, 
momentum and energy transferred through electron Green’s 
function are denoted. Line (2): Diagrammatic representation 
of Bethe-Salpeter equation for Cooperon Co-o-'; 77 r 7 '(s, q; w', tj). 
The parametrization of the Cooperon frequencies is shown to¬ 
gether with spin states. Line (3): Keldysh structure of the 
equation’s kernel. 


has been found to be shorter than the dephasing time. 
This permits us to study only valley-singlet Cooperons 
and ignore a possible sublattice and valley structure of 
elastic (spin-less) disorder in graphene, hence, to derive 
the form of the Cooperon equations affected by spin- 
flip scattering using a simplified model for electrons with 
structure-less (on the unit cell scale) Bloch functions. 

The diagrammatic form of the Bethe-Salpeter equation 
for the Cooperon matrix C is shown in Fig. 1(2). Thick 
solid lines in Fig. 1 denote disorder-averaged single¬ 
electron Green’s functions Ga (obtained from the solu¬ 
tion of Dyson equation in Fig. 2), 

= {ea-VF£.p+gegBB\\a +T-^))~^ ■, 

= {I-2np{e,)){G^-G^y, 

where 

T-l = [1 - 2crA/i(l - 2nF{£-a))] T^Y 
Afn = (5r)/5(5 + l); ^p^vf{\p\-Pf), 
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FIG. 2. Line (4): Diagrammatic representation of Dyson 
equation for single electron Green’s function. The last dia¬ 
gram contains the selfconsistent self-energy calculated in the 
main order of IjvppT expansion. Line (5): Keldysh structure 
of the self-energy. 


and <7 is electron’s spin projection on the direction of 
in-plane magnetic field. Here 


n^U HjJSa HjJS a* njJS*<J 

(7,E^+C0 )],E^ Cr,E^ (7,E^ i,E^+CO '['.Ej IAj 


( 6 ) 
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c) d) 


(5r) = \Z{a)\-^dlZ{a), 
where a = gipLBB]i^/kT, and 

Z{a)= Y. 

s,=-s 

is a statistical sum for the spin ensemble of paramagnetic 
scatterers. 

The disorder correlation function is represented in 
Figs. 1 and 2 by the dashed lines. Here, we assume short- 
range spatial correlations and include the following ele¬ 
ments in the description of disorder listed in Fig. 3: 

(a) spatial correlator has only Keldysh component 

{2'k^t)~^5{uj)5{y — r') 

for the short range spinless disorder (the first term in 
H(r)), which is parametrised by the momentum relax¬ 
ation rate, = 2'K^nuU‘^ and the electron density of 
states 7 ; 

(b) correlators of z-spin components for paramagnetic de¬ 
fects that characterise spin-dependent scattering of elec¬ 
tron without spin flip. 


FIG. 3. Line (6): The Keldysh element of the diagram 
corresponding to the two-electrons scattering on the impu¬ 
rity (after averaging over disorder).: (a) dominating spin- 
independent scattering, (b) elastic spin-dependent without 
spin-flip, and (c,d) spin-flip scattering. Line (7): The re¬ 
tarded (advanced) element of the diagram corresponding to 
the two-electrons spin-flip scattering on the impurity (c,d). 
The parameters of incoming and outgoing Green’s function 
are denoted. 

where correlator D{t,t') = ^KS+{t)S-{t')), written as 
a matrix in Keldysh space [2|, has the following compo¬ 
nents 0, 

= 2i{S^){uj - g,pBB ± f0)"\ 

D^{u:) = 47r[5(5 + 1) - (52)]d(a; - 

The implementation of Keldysh functions of the de¬ 
fects’ spins in the perturbation theory diagrams is shown 
in Figs. 1 and 2. Summing all three possible combina¬ 
tions of Keldysh components noted in Fig. 1, 

G^D^G^, and G^D^G^), where the frequency argu¬ 
ment of the spin corellator is -b 

giPBB), we get the contribution to the kernel of Bethe- 
Salpeter equation for the Cooperon, including all spin 
exchange processes (Fig. 3 c) and d)): 


T~^ = 2'iT^njj'^, which is independent on the positions 
t and t' on the Kendysh contour and therefore also has 
only Keldysh component^ Fourier transform of which is 


27r7Ts 


(5(a;)d(r 


(c,d) spin-flip correlators for paramagnetic defects, de¬ 
fined on the Keldysh time contour. 


-^(T^5+(t)5_(0)<5(r-r'), 

L'K'^Ts 


^+ [1 _ 2nF{e. + u:)]{D^ - D^)) 

= +l)-{Sl)-[l-2nF{e.)]{sS)K^). 

We note that in the approximation where only the 
Cooperons C(e, q; w', w) oc d(w') are taken into account, 
the energy transferred through the impurity spin can be 
only where the sign depends on whether the 

impurity transfers spin 1 or —1. 

After this, the equation for the Cooperon takes the 
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form 


(^-ihdr - C(e, r) = ^^Sir); 

R ( 7 ( 7 ' \ r ] r ]' — ^ ar )^ a ' r }' ^ a ' )/^ ^ 

+ 5a+a',r)+-n' (-4(TCr'r7^A /2 + |cr - Jyka ^ + 

+ i(cr-cr')(5e-5i)MB^), (6) 


where D = v^t/2 is a diffusion coefficient, C(e,r) is 
the Fourier transform of C(e, q) and = (0,i?_La;) is 
the vector potential of the perpendicular magnetic field, 
which is assumed to be small {B± B\^). Note that here 

we have already used a simplification of the Cooperon 
equations in graphene by assuming that both decoher¬ 
ence time and spin flip time in the material are longer 
than the intervalley scattering time, t^, ^ Tj„, so that 

these Cooperons have to be understood as valley-singlet 
Cooperons 0 (one does not need to involve valley-triplet 
Cooperons in the analysis). 

Diagonalization of the matrix R produces Cooperons 
decoupled from all others with decay rates 

= [1 - M 2 T Mi(l - 2nj^(£:p))] T.-i + (7) 

where takes into account that spin-flip scattering re¬ 
quires energy transfer to the scatterer. At the same time, 
spin precession as well as spin-flip processes mix up two 
Cooperons, Cq.o and Ci^o- This mixing generates com¬ 
bined modes, which relax with the decoherence rates 

= [1 + M2 - Mi{nF{e+) - nF(e-))] t~^ 

- {9e - g^y^llBj/^f + ^ ^ 

Note that, at B\\ = 0, (Sz) = 0 and (S^) = ^S{S + 1), so 

that ~ '^cT- ~ ''’ 0 "+ = 

corresponding to the relaxation rates of triplet and sin¬ 
glet Cooperons in the spin-memory regime. At B^\ ^ Bt, 
reflecting the restoration of phase coherence 
of diffusion of spin-polarized electrons. 

Together, the four Cooperon modes obtained from 
Eq. (O and substituted in Eqs. (I2|) yield the WL cor¬ 


rection to the conductivity at B^ = 0 




a—± 


In -|- Aa In — 

Tiij 




(9) 


1/A± = (Tp 1 -T. 


-1 
T I 


where Tiv is the intervalley scattering time. 

The magnetoconductance (defined as conductance 
variation as a function of B± for fixed B^ takes the form 


t{B±,Bi\) -a{0,Bi\) = - 


2Trh 


de ^ n'p{ea) 


a—± 


X 




F{z) = In^; -I- 1 /) 




. ( 10 ) 


71/3,a 


h/Ae 
Drp^a. ’ 


where '0 is digamma function. 

By using the expansion F{z) k, z^I2A -|- 0[z^], we 
obtain the formula for the magnetoresistance curvature 
(which gives a direct measure for the decoherence rate), 


j f DjXaX^ -I- Tq- Dtq,+ Dtq^^ \ 

^ \ V ^/e / ^/e fi/e j 


( 11 ) 


This result was used in the main text to reproduce the ex¬ 
perimentally observed non-monotonic dependence of MR 
curvature on the in-plane magentic field and to explain 
the anomaly in the ’decoherence rate’ that take place 
when l^e - -- hr-^. 
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